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We introduce a random matrix model with the symmetries of QCD at finite temperature and 
chemical potentials for baryon number and isospin. We analyze the phase diagram of this model 
in the chemical potential plane for different temperatures and quark masses. We find a rich phase 
structure with five different phases separated by both first and second order lines. The phases 
are characterized by the pion condensate and the chiral condensate for each of the flavors. In 
agreement with lattice simulations, we find that in the phase with zero pion condensate the critical 
temperature depends in the same way on the baryon number chemical potential and on the isospin 
chemical potential. At nonzero quark mass, we remarkably find that the critical end point at nonzero 
temperature and baryon chemical potential is split in two by an arbitrarily small isospin chemical 
potential. As a consequence, there are two crossovers that separate the hadronic phase from the 
quark-gluon plasma phase at high temperature. Detailed analytical results are obtained at zero 
temperature and in the chiral limit. 
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I. INTRODUCTION 



Currently, there is a strong interest in exploring the 
phase diagram of QCD at finite baryon density. A large 
number of possible phases have been suggested for QCD 
at finite baryon density (see 0,0 for a review). However, 
at this moment, the existence of any of these phases has 
been confirmed neither by first principle calculations nor 
by the phenomenology of heavy ion collisions and neutron 
stars. Because of the phase of the fermion determinant, 
standard Monte-Carlo simulations are onlypossible for 
small values of the chemical potential [E EL IE IE IE IE E| • 
Neutron stars are probably the most likely candidates 
for high baryon density physics, but they are hard to 
observe and only few parameters can be measured accu- 
rately. Relativistic heavy ion collisions explore the region 
of low baryon density and high temperature but give a 
complex picture of QCD at finite density. However, an 
experimental observation of a tricritical point might be 
within the realm of possibilities. Such a tricritical point 
was predicted on the basis of effective potentials llflj and 
simplified models such as random matrix models [ljj and 
Nambu-type models |12(. Furthermore, in both neutron 
stars and relativistic heavy ion collisions, the isospin den- 
sity is different from zero. It is therefore important both 
phenomenologically and theoretically to study the influ- 
ence of isospin on the phase diagram of QCD at nonzero 
baryon chemical potential. 

Our main goal is to study the phase diagram of 
QCD at nonzero temperature and chemical potentials 
for baryon number and for isospin. The phase dia- 
gram of QCD with any of these external parameters 
equal to zero has already been studied in a variety of 
ways. In particular, because lattice simulations are pos- 
sible at nonzero isospin chemical potential |l3l Il4j , the 
plane of zero baryon chemical potential has been un- 
derstood best. At zero temperature, we expect a sec- 



ond order phase transition to a phase of condensed pi- 
ons at an isospin chemical potential equal to half the 
pion mass. Such type of transition occurs in any QCD- 
like theory with a chemical potential for the c harg e of 
a Goldstone boson [IE EE llS O, O IS3, ISH H^l 
EE Ell • This prediction from effective Lagrangians 
has been confirmed by numerous lattice QCD simulations 
El ElEEEl EE El El EE E| From lattice simula- 
tions |l4l|3j| and a one-loop analysis of the effective La- 
grangian EE EH it also follows that the second order line 
changes into a first order line at a tricritical point with 
critical chemical potential and temperature on the scale 
of the pion mass. At nonzero baryon chemical potential, 
analytical results have been obtained for asymptotically 
large values of the chemical potential where QCD can 
be analyzed perturbatively EE Ei| • Lattice QCD simu- 
lations are only reliable for small values of the chemical 
potential [E EI @ • This leaves us with the bulk of the 
/iB-T-plane which could only be analyzed in simplified 
models such as Nambu Jona-Lasinio models |l'2l l37| , in- 
stanton liquid models |38| and random matrix models 

El 

Random matrix models were introduced in the con- 
text of QCD to describe the correlations of the low-lying 
eigenvalues of the QCD Dirac operator [39l |4fJ . It was 
shown that these models are equivalent to the mass term 
of a suitably chosen chiral Lagrangian which is deter- 
mined uniquely by the s ymm etries of the underlying mi- 
croscopic theory 0,EEME0|- Therefore, in the chiral 
limit, chiral Random Matrix Theories provide an exact 
analytical description of the low-lying Dirac spectrum. 

In this article we use a Random Matrix Theory as a 
schematic model for phase transitions in QCD at nonzero 
temperature and chemical potentials. Such a model was 
first introduced in |45j| to describe the chiral phase tran- 
sition in QCD at nonzero temperature. Even more suc- 
cessful was the application of Random Matrix Theory to 
QCD at nonzero baryon chemical potential. First, the 
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failure of the quenched approximation was explained an- 
alytically ^(|. Second, a tricritical point was found in a 
model at nonzero baryon chemical potential and temper- 
ature 11]. Third, algorithms for QCD at finite density 
could be investigated in detail |47ll4q . Fourth, the static 
part of effective Lagrangians for QCD with a chemical po- 
tential for the charge of Goldstone bosons can be derived 
from a chiral random matrix model |17| . In spite of the 
schematic nature of the random matrix model, we hope 
that it will teach us more about the plethora of possible 
phases that may occur in QCD. 

The organization of this paper is as follows. In sec- 
tion [n] we deduce from general arguments the expected 
features of the QCD phase diagram for nonzero temper- 
ature, baryon and isospin chemical potentials. In sec- 
tion IIIII we introduce our random matrix model. In sec- 
tion llVI we derive an effective partition function in terms 
of the meson fields and reproduce the mean field results 
obtained from a chiral Lagrangian. The phase diagram 
resulting from the random matrix model is obtained in 
section Ivl Concluding remarks are made in section IVll 



II. QCD AT NONZERO CHEMICAL 
POTENTIALS AND TEMPERATURE 

The QCD partition function at nonzero temperature 
and a chemical potential for each quark flavor is given by 

N f 

(JJdet(£> + m/+M/7o)) (2.1) 
/=i 

where the Euclidean Dirac operator is given by D = 
liii^n + i-^-ti) with 7 M the Euclidean 7-matrices, and 
is an SU(N C ) valued gauge potential. The quark masses 
are denoted by rrif, and pf is the chemical potential for 
each flavor. The average is over the Euclidean Yang- 
Mills action. Below we mainly focus on QCD with two 
flavors and nonzero baryon number and isospin chemical 
potential. In that case, the baryon number and isospin 
chemical potential are defined by 

Ms = 2 (Ml + M2), (2.2) 

MJ = \(Pi - M2)- (2.3) 

Before discussing the possible phases of random matrix 
models with the symmetries of the QCD partition func- 
tion, we first make some general remarks on its phase di- 
agram at nonzero temperature, isospin and baryon chem- 
ical potentials. 

The case of pb =0 and m n , T <§C Aqcd can be de- 
scribed in terms of a chiral Lagrangian. This Lagrangian 
has been analyzed to one-loop order [2^, [25| . At low T, 
a second order phase transition to a pion condensation 
phase was found at pi — where m w is the physical 

pion mass. For pi > m n /2 the chiral condensates rotate 



into a pion condensate, denoted by p, and approach zero 
for fij 3> rriir/2. Since a baryon chemical potential only 
excites states for p b larger than the nucleon mass, we ex- 
pect this second order line to persist in the pi-ps -plane 
and to be parallel to the pb axis. 

The pj — plane was discussed in detail in 11]. We 
expect a region with broken chiral symmetry (with chiral 
condensates (uu) and (dd) both nonzero) separated from 
the region of unbroken chiral symmetry by a first order 
curve, denoted by p c (T), from the tricritical point to the 
T = axis and a second order curve from the tricritical 
point to the ps = axis. 

For pj and \ib both nonzero, there are eight possible 
phases with either of the chiral condensates, (uu), or (dd), 
or the pion condensate p — ■^((u'jsd) — (chfsu)) equal to 
zero or not. Since the chemical potentials for the two 
flavors are different, there is no reason to expect that 
(uu) = (dd). 

In the limit p\ S> Aqcd, one flavor decouples and we 
are in a situation with only one flavor at nonzero chemical 
potential. In this case, we expect that (uu) = and 
(dd) ^ for pi < p c (T) but vanishes across the first 
order transition curve for P2 > l^c(T). 

For pb = 0, we have that pi = —pi- Using that 

det(£> + to - 1U170) = det*(D + m + pu ), (2.4) 

we find that for equal quark masses the partition function 
(|2.1() is the phase quenched partition function for two 
flavors [l| . 

III. RANDOM MATRIX MODEL 

In this article we study a random matrix model for 
QCD at nonzero chemical potentials and temperature. 
The idea is to replace the matrix elements of the Dirac 
operator by Gaussian random variables subject only to 
the global symmetries of the QCD partition function. 
The dependence on the temperature and chemical poten- 
tials enters through external fields structured according 
to these symmetries. 

Our guiding principle for constructing a random ma- 
trix model is dictated by the global symmetries of the 
QCD partition function. At zero temperature and chem- 
ical potential, this amounts to replacing the matrix el- 
ements of the Dirac operator by Gaussian random vari- 
ables subject to these global symmetries. The external 
field Q(pf,T) representing temperature and chemical po- 
tentials is introduced according to the following criteria: 

• The chemical potential breaks the global flavor 
symmetry in the same way as in the QCD parti- 
tion function. 

• The temperature field does not break global flavor 
symmetries. 

• For an anti-hermitian Dirac operator, the temper- 
ature field is anti-hermitian, whereas the chemical 
potential field is hermitian. 
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• The eigenvalues of the external field are fif ± inT. 
In this article we only consider the case n = 1. 

For two flavors, the Dirac operator of the random ma- 
trix model is given by 



/ mi A W + uj(T) + fi 1 \_ 

-A m 2 W + u(T)+(i2 

-Wl+uj(T) + m mi -A 

\ -W 1 + w(T) + fi 2 A m 2 / 

(3.1) 



We have also included a pion condensate source term, 
which in QCD is given by 



i\lpj5T2lp, 



(3.2) 



where the Pauli matrix r 2 acts in flavor space. The ma- 
trix elements of the n x n matrix W are complex with 
probability distribution given by 

P(W) = exp(-nG 2 TrVyVF t ). (3.3) 

The temperature field given by the matrix 



w(T) 



iT 
-iT 



(3.4) 



includes only the two lowest Matsubara frequencies. At 
zero chemical potentials, this temperature dependence 
will result in a second order phase transition along the 
temperature axis |45| . If we write the determinant as a 
Grassmann integral, the partition function of our model 
is given by 



Z = J VWY[d^ d^ S P(W)exp[-^2il) f 



Dip 



(3.5) 



where D is the Dirac matrix given in (|3~T|) . In the ther- 
modynamic limit the partition function is a function of 
mi, m 2 , A, hi, \jlb and T, but for brevity we will not 
display its arguments. 

Other ty pes of random matrix models have also been 
considered [4!3, |5(J, |5l|, |52| • However, none of these mod- 
els have been studied at nonzero isospin chemical poten- 
tial. We mention models with a random gauge potential. 
In that case, the matrix W has the spin and color struc- 
ture of the usual Dirac operator, 



W -> ia„T k A*, a v = (-i, <r k ) 



(3.6) 



but with Ajj a Gaussian nxn random matrix. Such types 
of Dirac operators have the same spectral properties |53| 
as the Dirac operator in l|3.5|l and lead to a similar phase 
diagram. 



IV. EFFECTIVE PARTITION FUNCTION 

Because of the unitary invariance of the random matrix 
models, the partition function can be rewritten in terms 
of invariant degrees of freedom only. Below we rewrite 
the partition functions introduced in Section fllll in terms 
of these effective degrees of freedom. 

We consider the random matrix model for QCD at 
nonzero temperature, baryon and isospin chemical po- 
tentials given in 13.5(1 . The Gaussian integration over the 
matrix elements of W can be performed trivially. The re- 
sulting four-fermion interaction is decoupled by means of 
a Hubbard-Stratonovich transformation at the expense of 
introducing mesonic degrees of freedom. After perform- 
ing the Grassmann integrations, the partition function 
(|3.5|) can thus be written as 



Z = J X>Aexp(-£(A, A 1 ")), 



(4.1) 



where 



C = uG 2 Tt(A - M f ){A^ - M) - -TrlogQ' 

= nG 2 Tr(A - M r ){A*< - M) ~ ^Tr log Q^Q, 

(4.2) 

and A is an arbitrary complex Nf x Nf matrix. The 
determinant of the 4iV/ x 47V/ matrix Q' , given by 

A iT + ^B +(J,iI 3 

A -iT + HB+Hih 

iT + n B +nih At 

-iT + HB+Vih At 

(4.3) 

with ^3 = diag(l, —1), factors into the determinant of 
A iT + fi B + Hih 



iT + hb+ Vih 



At 



(4.4) 



and its Hermitian conjugate. The mass matrix is given 
by 



M = 



mi —A 
A m 2 



(4.5) 



By shifting A, we have absorbed the dependence on M 
into the quadratic term. 

At zero temperature and chemical potentials, the chi- 
ral random matrix partition function is equivalent to the 
zero momentum part of the QCD chiral Lagrangian. We 
will now show that this is also the case for the chiral La- 
grangian that can be derived from (|4.1|) . To derive this 
result, we use the power counting scheme that is used in 
the construction of the chiral Lagrangian: /ij, ~ y/m 
and %/A are of the same order ^(| . We thus expand the 
Random Matrix Theory effective Lagrangian to first or- 
der in m and A and second order in fij about the saddle 
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point obtained for /ij = m = A = 0. The saddle point 
equation given by 

G 2 ((AA^ + T 2 - p 2 B f + Ap 2 B T 2 )A = (AA* + T 2 - p 2 B )A 

(4.6) 

has two solutions, 

A = 0, (4.7) 

or the solution of 
G 2 ((AA^ +T 2 - p 2 B ) 2 + 4/i|T 2 ) 



(AA^ +T 2 - n B ). 



(4.8) 

For hi = it is a natural assumption that the flavor sym- 
metry is not spontaneously broken. We can parameterize 
A in the broken phase as 



.4 



^cr(pi s ,T)E, 



where £ is a unitary matrix and 



1 1 



*(J*B,T)= [- + -^1-(AG^bT) 2 -T 2 G 2 



(4.9) 

p%g 2 

(4.10) 



V. PHASE DIAGRAM 

As was argued in 0, Hi| . the free energy (|4.12|) is 
completely determined by the transformation properties 
of the QCD partition function. Since the random matrix 
model has the same global transformation properties as 
the QCD partition function, we thus find the same low- 
energy limit. In this section we analyze the random ma- 
trix partition function beyond this universal domain. In 
our model we will be able to study the partition func- 
tion to all orders in the mass, chemical potential and 
temperature. We will show that the inclusion of such 
non-perturbative contributions alters the nature of the 
phase transition. 



A. Observables 

We consider three different observables, the chiral 
condensates (uu) and (dd), and the pion condensate 
i/2^ ((^75^) — (dj5 u ))- They can be expressed in terms of 
derivatives of the partition function, 



— d mi \ogZ 
2n 



cir 



is a solution of the saddle point equation l|4.8|l . Using 
the ansatz (|4.9|l . the inverse of the matrix Q for fj,j = 
is given by 



G 2 -<A* 1+ Au) 



Q- 



l 



a 2 /G 2 - {fiB + iTY 



At 
~Pb - iT 



A 



(dd) = 



(4.11) 



k 9 > 

G 2 



m 2 fog Z 
1 



cir 



(A* 22 + A 22 ) 



mi 



m 2 



(5.1) 



(5.2) 



Inserting this result and the parameterization (|4.9|l in 
the chiral expansion of the Lagrangian l|4.2|) . one easily 
derives 

C = c {n B) ni,T)-nGa(fi B ,T)Tr(M^ + M f £) 
+nriG 2 <j 2 (fi B ,T)c 2 {fi B ,T)Tii^I 3 ZI 3 ), 

(4.12) 

where cq(pb, T) is independent of £ and 

{AG 2 [i B T) 2 
Mai(p B ,T)-G2(p B -T2)y 

(4.13) 



-((wy 5 d) - (d'jBu)) 

= -U A logZ e 
4n 



-(A 12 +Al 2 -A 21 -A* 21 )-\).(5.3) 



C2(MB,T) = 1 - 



This effective Lagrangian coincides with the zero mo- 
mentum part of the leading-order Chiral Lagrangian 
at zero temperature and baryon chemical potential de- 
rived in [iT], [l8j based on the symmetries of QCD. A 
transition to a pion condensation phase takes place at 
p 2 c = m / (2Ga(/j,B,T)c 2 (^B,T)). The pion condensates 
vanishes for /x/ < pi >c - For pi > pi tC the chiral conden- 
sate rotates into a pion condensate but the sum of the 
squares of the chiral condensate and the pion conden- 
sate remains constant. The temperature and the baryon 
chemical potential affect both the magnitude and the ori- 
entation of the condensates. 



G 



The expectation values of the diagonal matrix elements 
of A can be interpreted as the chiral condensates, whereas 
its off-diagonal elements represent the pion condensate. 



B. Effective potential 

In the large-n limit, the partition function can be cal- 
culated by a saddle point approximation. To solve the 
saddle point equations, we make an ansatz for the matrix 
A. Since we have two independent chemical potentials, 
the chiral condensates are not necessarily equal, whereas 
for a sufficiently large isospin chemical potential we ex- 
pect a pion condensate. Based on the expressions for the 
chiral condensate and the pion condensate, we make the 
following ansatz for A 



A = 



0i p 

-P o 2 



(5.4) 
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In this parameterization, the different condensates are 
given by (uu) — G 2 (ui — mi), (dd) = G 2 (cr 2 — m 2 ), and 
^((uj^d) — (J75U)) = G 2 (p — A). Using this ansatz, we 
obtain the effective potential 

-C = G 2 (((Ti - m) 2 + (<r 2 - m) 2 + 2(p - A) 2 ) 
n 

"5 E l0 § [((^ + 0*1 ± iT ))( ff 2 - (M2 ±iT))+p 2 ) 

x ((en - ( M i ± iT))(a 2 + (p 2 ± tT)) + p 2 )] . (5.5) 

From here on, we set vn\ = m 2 = m. For p — this 
effective potential is a function of p 2 . In particular, this 
means that its dependence on pi at ps = is the same 
as its dependence on ps at pj = 0. This implies that for 
the critical temperature we have the relation [3, 0] 

rc(/*j)U=o = r c (MU=o fOT P = 0- (5-6) 

The fermion determinant of the theory with /is = and 
equal quark masses is equal to the fermion determinant 
of the phase quenched partition function. We thus ex- 
pect that the phase quenched approximation works in 
the phase where the pion condensate vanishes. 

To find the phase structure of our partition function, 
we have to solve the saddle point equations 

75 — = 0, — = 0, ^-=0. 5.7) 

0(Ji 0(J2 Op 

We will treat the following cases analytically: the chiral 
limit (m = 0) at zero T and at finite T, and the case 
m ^ at T — 0. We will calculate the phase diagram 
for T / and m ^ by numerically minimizing the 
effective potential. In addition, some analytical results 
are obtained for large quark mass m. 



C. Chiral limit at T = 

In the case of vanishing quark mass, zero diquark 
source, and zero temperature, the effective potential sim- 
plifies to 

-C = G 2 {a 2 + a 2 2 +2p 2 ) 

n 

1 • 2 

~2 lo g (( CT i + ^i)( CT 2 - M2) + P 2 ) 

--log ((<7! - p x ){a 2 + p 2 ) + P 2 ) 2 ■ (5.8) 

In the chiral limit, chiral symmetry is broken sponta- 
neously. As soon as the isospin chemical potential is 
switched on, the chiral condensate rotates into a pion 
condensate. Therefore, no phase exists where both con- 
densates are nonzero. The saddle point equation for p 
has two possible solutions: p = and p 7^ 0. We first 
consider the case p = 0. 



For p — 0, the effective potential separates into the 
sum of free energies for ci and a 2 , 

\C = G 2 o)- X -\o%{o)~p 2 f ) 2 . (5.9) 

/=1,2 

The saddle point equations given by 

°7 (^ 2 (°/ ~ M/) _ l) = 0, / = 1,2, (5.10) 

have the solutions 

a f = 0, / = 1,2, (5.11) 

o) = / = 1,2. (5.12) 

The contribution to the free energy from one flavor for 
these solutions is given by 

n f = -logp 2 , (5.13) 
n f = l + logG 2 +p}G 2 , (5.14) 

respectively. The full free energy is a sum over the con- 
tributions from both flavors. The two solutions are sep- 
arated by a first order phase transition line where their 
free energy is equal, 

l + ^G 2 + log(/4G 2 ) = 0. (5.15) 

The solution of this transcendental equation is given by 
PfG = p c G « 0.527697 0. At this point a first or- 
der phase transition to a phase with af — takes place. 
In the /ii-^2-plane, we can thus distinguish four differ- 
ent phases with nonzero condensates in strips along the 
chemical potential axes. The strips overlap in the cen- 
ter and form a region where both chiral condensates are 
nonzero. 

At zero quark mass, the Goldstone bosons are mass- 
less, and the critical value of the chemical potential for 
pion condensation is pi — 0. For pj > 0, the chiral 
condensates rotate completely into a pion condensate, so 
that a 1 =02 =0. The effective potential for p is given 
by 

-C = 2G 2 p 2 -\og(p 2 ~p 1 p 2 ) 2 . (5.16) 
11 

The saddle point equation given by 

p(G 2 ( P 2 -p 1 p 2 )-l) = (5.17) 
has again two solutions, 

p = 0, 

P 2 = -g2 +MiA*2. (5.18) 

A second order transition line is given by the hyperbola 
p 2 = in the quadrants where p\P2 < 0. The free energy 
of this phase is given by 

n p = 2(l + logG 2 + p 1 p 2 G 2 ). (5.19) 
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FIG. 1: Phase diagram of the RMT model for three colors in 
the chiral limit (m = 0) at zero temperature. Solid (dashed) 
lines are lines of first (second) order phase transitions. Except 
for the four corner regions, the different phases are marked by 
the nonvanishing condensate. In the four corner regions we 
have that p = o\ — 02 — 0. 



D. Chiral limit at T / 

At nonzero temperature, we analyze the phase struc- 
ture in the same way as we did for T = 0, and, initially, 
we also find the same phases. 

In a phase with vanishing pion condensate, the effective 
potential again separates into a sum over contributions 
of two different flavors, 

l -L = £ G 2 o) - \ log((a 2 + T 2 ) 2 + A, 2 T 2 ). 



/=1,2 



(5.20) 



This free energy was studied in for the case of zero 
isospin chemical potential. For each of the two flavors 
(/ = 1, 2), the saddle point equation 



fa 



2<y 2G 2 
T 2 



(4 



T 2 )a 2 



2x2 



G 2 



(»j+T 2 ) 



f 
= 



has the solutions 
a s = 0, 



(5.21) 



(5.22) 



-1 



1 

2G 2 



+ p) 



T 2 ± 



1 

2G 2 



l-{AG 2 n f Tf 



The free energy for each flavor of these solutions is equal 
to 



Finally, for /ii — /i2, the saddle point equations allow 
a solution with p ^ and ui = a 2 / 0. However, the free 
energy of this solution is higher than that of the solution 
with (Ti = (72 7^ and p — 0. Physically this is clear, 
since the pion condensate is expected to vanish for zero 
isospin chemical potential. 

Comparing the free energy of the pion condensation 
phase to those of the chiral condensation phases and the 
chiral restored phase, we can determine the remaining 
first order phase transition lines. The resulting phase 
diagram is shown in Fig. ^ We find a region of pion 
condensation in the center of the phase diagram. It is 
bounded by first order phase transitions towards phases 
with nonzero chiral condensate for one flavor, and by 
second order transition lines towards the chiral restored 
phase. The chiral condensation phases form the arms 
of a cross along the chemical potential axes. Since a\ 
is independent of /12, it is nonzero in a strip along the 
/i2-axis, and the same applies for the other flavor. The 
first order lines intersect in the two points (piG,p2G) w 
(±0.527697, ±0.527697). The intersection points of the 
second order lines with the first order transitions be- 
tween the pion and chiral condensation phases are at 
(/iiG,/i 2 G) w (±0.527697, Tl-895025), and at the two 
points obtained by interchanging the values for p,\ and 

M2- 



log(Ai 2 



/ 

logG 2 



(5.23) 



G 2 (pj-T 2 )±^Jl 



1 / 1 1 
--log - ± - 

2 S V 2 2 



1 - (4G 2 p f ry 



(4G 2 /i/T) 2 
(5.24) 



respectively. For (4G 2 /i/T) 2 < 1, the solution with the 
negative branch of the square root will be discarded, since 
it does not minimize the free energy. The second order 
phase transition line is given by the condition that the 
two solutions coincide. One easily derives 



{A 



T 



2\2 



tf - -l j T G (fj, f 

Since there is always the solution er/ 







(5.25) 
0, there can be 



a first order transition when the coefficient of in the 
saddle point equation l|5.21l) becomes negative. A tricrit- 
ical point occurs where both the coefficient of <jj and 
in H5.21f> vanish. This results in the equations 

1 



/4 



2G 2 

rp2 1 /~l2 



G\p) 



■ T 



2 -,2 



0. 



(5.26) 



with solution given by 



2 G 2 - 



2,-<2 



T 2 G 



V2-1 
1 : 
V2 + 1 



/ = 1,2, 



(5.27) 
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Numerically, T 3 G « 0.776887, fi f>3 G « 0.321797, which 
was also obtained in 

In Fig. |21 we show the phase diagram for fixed temper- 
ature in the p \-p 2 chemical potential plane. The transi- 
tion lines for each flavor are straight and constant in the 
chemical potential for the other flavor. In this plane, the 
phase transition lines change in their entirety from first 
to second order when we pass the tricritical temperature 
(|5.27() from below. The critical chemical potential of the 
transition becomes smaller with increasing temperature. 
A baryon chemical potential that is large enough will de- 
stroy the pion condensate. 

Since the two chiral condensates are independent of one 
another and depend only on the chemical potential for 
the respective flavor, we again have four phases. A phase 
where chiral symmetry is restored (a± — <y% — 0), and 
phases where either one or both of the chiral condensates 
are nonzero. The free energies are simply given by the 
sum of the one-flavor free energies. 

At nonzero isospin chemical potential, we expect that 
in the limit of massless quarks the chiral condensates are 
completely rotated into a pion condensate. The effective 
potential for vanishing Of and nonzero pion condensate 
p becomes 

-C = 2G 2 p 2 -\og((p 2 -p 1 p 2 +T 2 ) 2 +T 2 (p 1 +p 2 ) 2 ). 
n 

(5.28) 

The saddle point equation reads 

G 2 P (p A - 2(^j + P!P2 ~ T 2 )p 2 + ^( M l/i2 - T 2 ) 

+ (^2 - T 2 ) 2 + T 2 {p l + p 2 ) 2 ) = 0. (5.29) 
It has the solutions 

p = 0, (5.30) 
P 2 = ^2 + - T 2 ± ^Lyfi - (2G2( Ml + p 2 )Tf. 

Again the solution with the negative branch of the square 
root has a larger free energy. A line of second order phase 
transitions is determined by the condition that the two 
solutions coincide: 

{piP2 - T 2 ) + G 2 (pip 2 - T 2 ) 2 + G 2 T 2 (p 1 + p 2 ) 2 = 0. 

(5.31) 

A first order phase transition may occur when the co- 
efficient of p 3 in (|5.29|) vanishes. However, we will see 
below that this happens in a region where solutions with 
nonzero chiral condensate have a lower free energy. 

In Fig.|3we show the phase diagram in the pi- p 2 -p\a,ne 
for zero quark mass and temperatures equal to TG = 0.3, 
TG = 0.5, TG = 0.6 and TG = 0.8. The first order 
lines in the phase diagram are obtained by combining 
the results for the free energies of the phases discussed 
above. 



The phase diagram at T = has been described in 
the previous section. With increasing temperature, the 
second order transition between the pion condensation 
phase and the chiral restored phase moves towards the 
origin. The effect of the temperature term on the phase 
diagram is a shrinking of the condensate phases. The 
critical chemical potential for the chiral restoration tran- 
sition decreases with temperature, so that the transition 
lines move toward the axes as well. 

At lower temperatures, the phase where both chiral 
condensates are nonzero simultaneously is always higher 
in energy than the pion condensation phase, and conse- 
quently it is not realized. At a temperature of TG rs 
0.548047, a first order transition between the pion con- 
densation phase and the phase with <j\ ^ 0, a 2 ^ 0, and 
p = emerges, and it appears at the intersection points 
of the first order transitions, p\G — p 2 G w 0.413485, 
around the line p\ = p 2 . The upper boundaries of this 
phase are always the transition lines where either of the 
two chiral condensates vanish. 

The position of the tricritical point in the pf-T-plane 
for any of the two chiral condensates cr/ is unaffected by 
the presence of the second chemical potential. There- 
fore, we find that at a temperature T 3 G — | y \[2 + 1 w 
0.776887, the phase transition lines between the chiral 
condensed phases and the chiral restored phase in the 
chemical potential plane become second order transition 
lines in their entirety. 

Above this temperature, the regions with nonzero chi- 
ral condensates along the p\- and p 2 -oxe& are bounded by 
second order transition lines to the chiral restored phase. 
We can observe the region where both phases with chi- 
ral condensation overlap. The pion condensation phase 
in the center is separated from the chiral condensation 
phases by first order transition lines, and from the chiral 
restored phase by a second order transition line. 

All condensation phases vanish at the same tempera- 
ture TG = 1. 



E. Zero temperature limit at nonzero quark mass 

Away from the chiral limit, we can only solve the sad- 
dle point equation analytically for T = 0. The results 
of this particular case will be discussed in this section. 
The phase diagram at nonzero quark mass is qualita- 
tively different from the phase diagram in the chiral limit. 
First, the chiral condensates, o\ and a 2 , are not good or- 
der parameters anymore, and, second, we expect a phase 
transition to a phase with nonzero pion condensate for 
pi = mn/2. However, in this case the chiral condensates 
are non-vanishing in the phase with p 0. 

The effective potential is given by 

-C = G 2 ((a 1 -m) 2 + {<j 2 -m) 2 + 2p 2 ) 
n 

1 2 
--log ((<7i +pi)(a 2 -M2) +P 2 ) 
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FIG. 2: Phase diagram of the RMT model in the chiral limit for different temperatures. Solid (dashed) lines are lines of first 
(second) order transitions. The values for the temperature are given in the figures. Except for the four corner regions, the phases 
are marked by condensates that do not vanish in the chiral limit. In the four corner regions we have that p = a\ — 02 = 0. For 
TG — 0.6, the phase with o\ ^ 0, 02 7^ has emerged. For TG = 0.8, the chiral restoration transitions have become second 
order transitions. The region in the center where both chiral condensation phases overlap and the chiral condensates for both 
flavors are nonzero can be seen clearly. All condensates vanish at TG = 1. 



1 log((a 1 -p 1 )(a 2 +p 2 ) + p 2 ) 2 . (5.32) 2G 2 (a 2 -m)- 



2 V J (Ol +/Xl)(<7 2 - + P 

(71 - Mi 

(0-1 - pi)(a 2 + p 2 ) + p 2 ' 



Obviously, it is symmetric under a simultaneous inter- +7 tt : 5-, (5.34) 



change of the chiral condensates a\ and a 2 and the two 

chemical potentials p\ and p 2 . The saddle point equa- 2G 2 p = — — 

tions are given by ' ' (01 + Mi)(o"2 - M2) + P 

P 



or 2 („ ml - °2 (<7i - pi)(a 2 + p 2 ) + p- 

{ai-m) — j — — — j 

(cti + pi)(a 2 - p 2 ) + p z 



,2 

(5.35) 



(7 2 + p 2 (533) As we see from the last equation in (I5.35|) . there is always 



(<7i — p\){(7 2 + p 2 ) + p 2 a solution with p — 0. We expect that this is the actual 
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minimum of the free energy below the critical chemical 
isospin potential for pion condensation. In this 
for m = 0, the saddle point equations for the two chiral 
condensates decouple and are given by 



G 2 (a f -m)(a 2 f -p? f )-a f = 0. 



(5.36) 



Although this third-order equation can be solved analyt- 
ically, it is more instructive to expand it in powers of to. 
To first order in to, we find for aj, / = 1, 2, 




The free energy of these solutions is given by 

n f = 1 + logG 2 + p)G 2 - 2mG (5.39) 



f 



+0(m 2 G 2 ), 
-\og(p 2 )-0(m 2 G 2 ), 



(5.40) 



respectively. The solutions l|5.37|l minimize the free en- 
ergy for small values of the chemical potential, and the 
solution (|5.38() minimizes the free energy for large values 
of the chemical potential. As in the chiral limit, we once 
again have four phases, where either chiral symmetry is 
broken spontaneously and the chiral order parameters 
G<jf are of 0(1), or where it is broken only explicitly 
and the chiral order parameters are of the order of the 
quark mass, 0(m). The free energies in these cases are 
given by the sum of the free energies for each flavor. In 
contrast to the case of the chiral limit to — 0, the phase 
with o\ 7^ 0, a 2 7^ appears at the center of the phase 
diagram. 

By matching the free energies for the phases with large 
and small values of the chiral condensates, we obtain the 
correction to the critical chemical potential due to the fi- 
nite quark mass m. The critical chemical potential shifts 
to 

m' /iC G = p c G (l + i + "^ G2 + 0(to 2 G 2 )) ,(5.41) 

where p c G is the result for the chiral limit, obtained in 
(ETTBT) . 

For to 0, both the pion condensate as well as the 
chiral condensates are nonzero in the phase with p ^ 0. 
We have to solve the full system of three saddle point 
equations. In fact, in this case the analytical solution is 
relatively simple. The two chiral condensates are related 
by the equation 



<J\ — cr 2 



Mi + M2 
i . 

Mi - M2 



(5.42) 



M 2 G 



1 

/ 




■ 

mG = 0.1 
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/ 
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2 -1 
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FIG. 3: Phase diagram of the RMT model for three colors 
at a value of the quark mass of mG = 0.1. Solid (dashed) 
lines are lines of first (second) order transitions. Except for 
the four corner regions, the different phases are marked by 
the condensate that does not vanish in the chiral limit. The 
condensates that are not displayed are of 0(ra). In the four 
corner regions we have that p = 0, cri = 0(m) and 02 = 
0(m). 



The solution for at, f = 1,2, is given by 



a I = m M/- 



+ 



2to 



1 



(Mi - M2) 2 ' G 2 {pi - ,U2) 2 - 4to 2 
The pion condensate then follows from 



(5.43) 



(<71 - to)(ct 2 - TO) + p 2 = ^2+MlM2- l 



TO *pip 2 



(Ml 



M2) 2 ' 
(5.44) 



The free energy of this solution is 



n(m,p 1 ,p 2 ) 



2 (1 + log G 2 + p±p 2 G 2 
^2 



TO G 



-to 2 G 2(mi 



M2J 



-- log 
2 B 



(Mi - M2) 
(Mi " 



M2J 



(Mi - M2) 2 - 4to 2 



.(5.45) 



For nonzero quark mass m, the phase in which p =/= 
does not extend to zero isospin chemical potential. For 
Pb = 0, the onset chemical potential follows by equating 
the free energy + of (|5.39l) to the free energy of the 
pion condensed phase given in eq. I|5.45|l . It is given by 
Mz c = 2% + C( m2 ) which identifies y / 2ro/G as the pion 
mass for pb = T = 0. 

By putting p 2 = in eq. I|5.44|l . we obtain the com- 
plete second order transition line that bounds the pion 
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condensation phase at low as well as high isospin chem- 
ical potential. Parametrized in terms of the baryon and 
isospin chemical potentials, it is given by 



i 



m 



1 



2G 2 p* - m 2 

o 

m 



,2^JJ 

2 

(J>I 



•Mb 



G 2 m| 



4G 4 (^ 2 - m 2 ) 2 



(5.46) 



For /is = we again find a critical chemical potential 
given by fj,j iC — rajlG + 0(m ). The phase diagram in 
the /ii-/X2-plane for mG = 0.1 and zero temperature is 
shown in Fig. 03 The qualitative difference to the mass- 
less case is the appearance of a region where both chiral 
condensates are nonzero in the center of the phase dia- 
gram. The dashed lines that border this region cross the 
/i/-axis (fis = 0) at fij — ±m w /2 and are roughly con- 
stant in fiB- They coincide in the chiral limit, and the 
central region becomes a phase of nonzero pion conden- 
sate and zero chiral condensates (see Fig. 0. For small 
quark masses m, the phase diagram at large values of 
the chemical potentials is almost unchanged in compar- 
ison to the phase diagram we found in the chiral limit. 
The transition lines are only shifted by small corrections 
of the order of the mass mG. 



F. Phase diagram at nonzero quark mass and 
temperature 

Finally, we study the phase diagram at nonzero quark 
mass and temperature. In this case it is no longer pos- 
sible to obtain an analytical solution of the saddle point 
equations. Instead, we determine the minimum of the 
free energy numerically. The results for the phases in the 
jLti-/i2-plarie for temperatures of TG = 0.3, TG = 0.5, 
TG = 0.6, TG = 0.8 are shown in Fig. H The quark 
mass used here is mG = 0.1. Below the critical tempera- 
ture, the phases are the same as at zero temperature (see 
Fig. |3J) and differ qualitatively from the chiral limit only 
in the central region (see Fig. 01. 

The phase diagram in the /xs-T-plane for \ij = has 
been studied in [lj . In the chiral limit, the chiral restora- 
tion transition extends as a second order line from the 
fiB = axis, changes order at a tricritical point, and in- 
tersects the T = axis as a line of first order transitions. 
For nonzero quark mass, the first order transition ends in 
a critical point, and the second order transition becomes 
a crossover. 

Fig. [3] shows the phase diagram in the /^s-T-plane at 
finite quark mass mG = 0.1 for zero isospin chemical 
potential and for fijG = 0.1. We observe that the first 
order curve splits into two first order curves that are sep- 
arted by 2fijG. This can be understood as follows. Below 
the threshold for pion condensation, the free energy sep- 
arates into a sum over the two flavors. For /i/ = 0, the 
chiral phase transition lines for both flavors coincide. A 
finite isospin chemical potential breaks the flavor symme- 



try, and the first order transition lines for the two flavors 
split and shift according to 



(5.47) 



where fi c (T) describes the transition line at /i/ =0. The 
critical temperature is not affected by the isospin chem- 
ical potential. 

We have seen in eq. (I5.6|l that in the phase with zero 
pion condensate the dependence of the critical temper- 
ature on the isospin chemical potential at zero baryon 
chemical potential is the same as the dependence of the 
critical temperature on the baryon chemical potential at 
zero isospin chemical potential. This suggests the possi- 
bility that for large values of the pion mass a tricritical 
point may appear in the /ib = plane. As we will see be- 
low, this turns out not to be the case. We first determine 
the domain in the /ib = plane where pion condensation 
occurs. 

Performing a similar analysis as below eq. I|5.41|l for 
the free energy at /ib =0 but nonzero m, T and /x/, 
we find that the region of nonzero pion condensate is 
bounded by the curve 

M - m 2 )G 2 - -m 2 - (/if + T 2 )(p 2 - m 2 ) 2 G 4 = 0. 

(5.48) 

For asymptotically large values of the quark mass, this 
curve reduces to the two expressions, 



(5.49) 



However, for \xb = the saddle point equations are 
also solved by p — 0. Then a first order transition takes 
place between the solutions with large mass expansion 
given by 



a, =a 2 =m±i^i-G 2 T 2 + 0(-^). (5.50) 

At the tricritical point these solutions merge with the 
extremum between them. From this condition we find 
that the position of the tricritical point is given by 



2,-<2 



TiG 



1 



1 



0(- 



1 



2 16m 2 G 2 m 3 G 3 ' 
p%, 3 G 2 = m 2 G 2 + 1- -\— + 0(- 1 



(5.51) 



2 8m 2 G 2 ' ~ v to 3 G 3/ ' 
and the value of a at the tricritical point is equal to 



1 



1 



^4^ + °W- (5 ' 52) 

One easily verifies from (|5.49|) that up to order l/m 2 G 2 
the tricritical point is inside the pion condensation re- 
gion, where the phase with a nonzero pion condensate is 
favored. Numerically, one finds that this is also the case 
for quark masses that are not asymptotically large. 
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FIG. 4: Phase diagram of the RMT model for three colors for a value of the quark mass of mG = 0.1. Values for the temperature 
are given in the figures. Solid (dashed) lines are lines of first (second) order transitions. Above the critical temperature, the 
chiral restoration transition becomes a crossover, denoted by dotted lines. Except for the four corner regions, the different 
phases are marked by the condensate that does not vanish in the chiral limit. The chiral condensates are not displayed when 
they are of 0(m). In the four corner regions we have that p = 0, o\ = (D(m) and ai = 0(m). 



VI. DISCUSSION 

Starting from a random matrix model at nonzero tem- 
perature and chemical potential for baryon number and 
isospin, we have obtained an effective potential for the 
matrix valued order parameter field. This order param- 
eter field arises naturally in this random matrix model 
that is based on the global symmetries of the QCD par- 
tition function. The expectation value of its diagonal ele- 
ments are the chiral condensates (uu) and (dd), whereas 
its off-diagonal elements give the pion condensate. To 
first order in and \x\ and zero baryon chemical poten- 



tial and temperature, the random matrix model coincides 
with the zero momentum part of the chiral Lagrangian 
that has been derived from QCD. However, the tricritical 
point found in lattice simulations |34| and in the chiral 
Lagrangian at nonzero temperature and isospin chemi- 
cal potential [2S| is not present in the random matrix 
model. We therefore conclude that the pion dynamics 
are important for the emergence of this tricritical point, 
as was suggested in j34(. 

Based on the effective potential, we have obtained a 
phase diagram for QCD at nonzero temperature, baryon 
and isospin chemical potentials. We have found a sur- 
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densates in our order parameter field. We find that close 
to the critical baryon chemical potential a small isospin 
chemical potential leads to a phase with (uu) of the order 
of Aq CD but (dd) reduced by a factor m/ Kqcd- 

In the phase with a vanishing pion condensate, the 
effective potential is an even function of the chemical po- 
tentials and separates into a sum of free energies for each 
of the two flavors. This has important consequences. 
Since the effective potential is even, the dependence of 
the partition function on /xg at /Ltj = is the same as 
its dependence on /xj at (1b = 0. Therefore, the phase 
diagram for baryon chemical potential smaller than the 
pion mass can be studied reliably by means of the phase 
quenched partition function. Because of the separability 
of the free energy, the critical curve for /xj = splits 
into two curves shifted by a distance of 2/xj. As illus- 
trated in Fig. 5, the structure of the phase diagram in 
the /x b-T '-plane is structurally altered by an arbitrarily 
small nonzero isospin chemical potential, even for mas- 
sive quarks. For a fixed /xj < we find that there 
are two first order phase transitions at small T when /xg 
is increased. Both first order lines end at the same tem- 
perature in critical endpoints with a separation propor- 
tional to /x/. The existence of two first order phase tran- 
sition lines and two critical endpoints might have very 
important consequences for the numerous phenomeno- 
logical systems where both \Ib and /xj are nonzero, such 
as neutron stars or heavy ion collision experiments. Fur- 
thermore, it has been shown that relativistic heavy ion 
collisions experiments might be sensitive to the critical 
end point in the fi s-T-plane for /xj = 0, |§3 ■ Our 
analysis shows that an increase in /xj results in a crit- 
ical end point with a lower value for /xs, thus making 
it easier to reach via heavy ion collision experiments. 
Our analysis also implies that two crossovers separate the 
quark-gluon plasma and the hadronic phase at small but 
nonzero baryon and isospin chemical potentials. There- 
fore the transition between these two phases should ap- 
pear smoother at fij ^ than at /xj = 0. These results 
have important phenomenological consequences. It is es- 
sential to confirm them by means of lattice QCD simula- 
tions or within other models. 



FIG. 5: Phase diagram in the UB-T-pl&ne for quark mass 
mG — 0.1 and an isospin chemical as shown in the label of the 
figure. A first order chiral restoration phase transition takes 
place at the full line that ends in the critical end point. For 
nonzero isospin chemical potential (lower figure) this curve is 
shifted in opposite directions for the chiral restoration transi- 
tion of (uu) and (dd) . The condensates that are not displayed 
are of 0(m). The dotted curves depict the crossover behav- 
ior. The temperature of the critical end point is not affected 
by the isospin chemical potential. 



prisingly rich phase diagram characterized by the con- 
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